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Abstract 

The effect of repeatedly smearing SU{3) gauge configurations is investigated. Six gauge 
actions (Wilson, Symanzik, Iwasaki, DBW2, Beinlich-Karsch-Laermann, Langfeld; com- 
bined with a direct S'C/(3)-overrelaxation step) and three smearings (APE, HYP, EXP) 
are compared. The impact on large Wilson loops is monitored, confirming the signal- 
to-noise prediction by Lepage. The fat-link definition of the "naive" topological charge 
proves most useful on improved action ensembles. 

1 Introduction 

Smearing of gauge-links has been used on many occasions as a powerful tool to reduce the 
UV-fluct nations of a gauge configuration [H El El IH E] • In principle, such a UV-filtering may 
be attempted either in the Monte Carlo process itself (by complementing the action with con- 
tributions from 14/1x2+2x1 even larger Wilson loops) or in the observable (where replacing 
the original links by smeared ones is just the simplest way to define a filtered observable). 
Normally, the motivation for the first option is not to tame the noise in standard observables, 
but to reduce cut-off effects. Nonetheless, an obvious question is whether there is an interplay 
between the action used to generate a pure gauge ensemble and the effectiveness of a smearing 
recipe to decrease the UV-fiuctuations in a given observable. 

Some expectation values depend in a rather sensitive way on the details of the smearing 
procedure. As an example. Table Q contains expectation values of a few Wilson loops without 
smearing and after two smearing recipes have been applied. Obviously, the plaquette and larger 
Wilson loops take rather different values in these three cases. However, since from the difference 
of the rows one may compute the heavy-quark (HQ) force, and such a physical observable is 
supposed to be unaltered by the smearing, the induced change in the W{R,T) must take place 
in a coherent manner for all R, T. Obviously, not only the core recipe (e.g. APE, HYP, EXP - 
see below) but also the parameters used and the number of iterations will have an impact on 
the quality of the result. Basic questions that one may have at this point include: 
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Table 1: Wilson loops {W{R,T)) with R,T e {1,...,4} on a 12^ lattice with Wilson glue at 
j3 = 6.0 (left), after one HYP-step with a = (0.3, 0.6, 0.75) and SU (3) -projection (middle) or the 
same a without any projection (right). Statistical errors range from 1 to 10 in the last digit. 
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1. For a given recipe, how can one identify useful combinations of parameter values and 
iteration levels, and is there a specific signature of doing "too much" smearing ? 

2. Optimizing parameter and iteration level w.r.t. a certain criterion, will the result depend 
on the gauge action used, or will it be universal for a given lattice spacing ? 

3. Is there a physics reason why some smearing recipes and/or parameters prove more effi- 
cient in damping UV-fluct nations than other ? 

The plan of this article is to first list and review the actions studied. Then the smearing 
procedures are discussed, shifting technical details of the S'?7(3)-projection to the appendix. 
The key physics point is an explicit test in the heavy-heavy case of the Lepage argument which 
establishes a relationship between the decrease of the signal in large Wilson loops and the self 
energy of an infinitely heavy quark. In order to give a fair comparison of different actions 
approximately matched ensembles must be generated and I describe how the associated (3- 
values have been found with little CPU costs via the local force. With these at hand, I study 
of the effect previewed in Table ^ combining all actions and smearing recipes. The result is 
that for Wilson loops at fixed lattice spacing the effect of smearing is almost independent of 
the gauge action used. Finally, I point out that on sufficiently fine lattices smearing allows for 
a cheap definition of the field-theoretic topological charge, most notably on improved glue. 

2 Action details 

We consider gauge actions involving rectangular Wilson loops up to maximum side length 2. 
The action involving nothing but the loop VT^^^ is the one due to Wilson jB] 

S'-' = /3 Y: l-\RetTWi::\x). (1) 
Under a change of the link U^{x) the acceptance is determined by the "local action" 

^lof = -f Retr(sixi(x)tf/,(x)) (2) 
where fi is not summed over and the staple is the sum of all 3-link paths from x to x+fi 

sl'^'ix) = Yl M^) U,{x+P) U,{x+fiy . (3) 

For future reference we define the rescaled Wilson staple Sj^^^ = ^ Sjj^^ where the prefactor 
simply refiects the 6 contributions to the sum in (jH)). 

The class involving the loops ly^^^ and pyix2+2xi ^ |yix2_|_p^2xi j^^^g ^^le generic form 

^'^'+'^^ = /?E (l-8ci)(l-^Retriy^V^i(^))+£i(2-^Retriy^if+2><i(x)) . (4) 

Under a change of U^{x) the acceptance is determined by the change of (no sum over /x) 

Slf^'x' = Retr([(l-8c0^;^^(x) + £i5f^+2xi(a:)]t U,{x)) (5) 
where 5*1x2+2x1 _ ^ix2_|_^2xi generates all planar 5- link paths from x to x+fi, viz. 

= E f/.(x)t/,(a;+i))t/^(x+2i))f/,(x+/i+i>)tf/,(x+/i)t (6) 

ST\x) = E U,ix)U^ix + iy)U^ix+fx+u)U,ix+2ix)^U^{x+ix)^ 
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The class of actions involving the Wilson loops 14^ and W has the generic form 

^2><2 = /?E (l-16c2)(l-^Retrl^^ir(^))+£2(l-^Retrl^^2f (x)) . (8) 



Under a change of U^{x) the acceptance is determined by the change of (no sum over /i) 

Slf = -^Retr([(l-16c2)5;xi(x) + C2Sf ^(x)]tf/,(x)) (9) 

where the new staple S^^^ is designed to generate those 7-link paths from x to x + fi which, 
together with f/^(x)^, form a 2x2 square, viz. 

Sl^'^x) = U,{x)U,{x + v)U^{x + 2v)U^{x + pi+2v)U,{x + 2fi+v)^U,{x + 2fi)^U^{x+fi)^ 

+ Uf,{x-iiy Uu{x-il)Uu{x-fx+u)Uf,{x-fi+2u)Uf,{x+20)U^{x+fi+uy (10) 

Specific choices for ci in the W^^^+W^^'^^'^^^ class include (ci = ci unless stated otherwise) 

{—1/12 {ci = ci,ci = ci/ul) Symanzik (tree-level, tadpole) [Zj 
-0.331 Iwasaki (KG) [HI (11) 

-1.409 DBW2 (RG) P 

with uo = (|Retr(iy^^^))^/^, while the C2 in the W^^^ + W'^^'^ ansatz that I am aware of are 

_J —1/48 (c2 = C2, C2 = C2/mo) Beinlich-Karsch-Lacrmann ftree-level. tadpole) jlO] 
1 —0.538 Langfeld (HQ potential isotropy) [TT] . 

(12) 

From (jni EI) and (fTTH) it follows that 5^^^+^^^ receives 18 contributions while 5^^^ is built 
from 12 terms. Disregarding the one-third matrix product needed to form the trace in tT{Sj^U^) 
(no sum) and counting only the multiplications for the staple, one gets 6 ■ 2 = 12 operations 
with the Wilson action, while the Symanzik-Iwasaki-DBW2 actions lead to 6 ■ 2 + 18 ■ 4 = 84 
matrix multiplications and for the Beinlich-Karsch-Laermann-Langfeld actions that figure is 
6 • 2 + 12 ■ 6 = 84, too. From such a counting one expects that the new actions are roughly 
a factor 7 more expensive (per sweep) than the Wilson gauge action, but of course the hope 
is that this investment pays off by e.g. a better scaling behavior. In any case it is clear that 
for the new actions the costs for a single-link update is dominated by the staple. In such a 
situation a multihit Metropolis is almost competitive with a heatbath, since 1 or 10 hits bear 
almost the same costs and after 10 hits the link is essentially in equilibrium. 

A point worth mentioning is the observation that at standard /5- values a "direct" overrelax- 
ation step in the full 5'[/(3)-group proves remarkably efficient, if properly implemented. With 
the Wilson gauge action the current link is proposed to be replaced according to 

U,{x) U'^ix) = Psui3){Sl^\x)} f/,(x)t Psui3){Sl^\x)}, (13) 

with Sjj^^ defined below (jS)) and Psuis) specified in the appendix. Unlike in SU{2) this proposal 
must be subject to an accept/reject step. It turns out that at /3 = 6.0 the acceptance rate in 
this over relaxation step is 99%. However, if we would keep the proposal ()13|) with the improved 
actions, the acceptance rate would degrade severely. Fortunately, the generalized proposal 

U,ix) ^ U'^ix)=Psui3){S,ix)}U,ix)^Psui3){S,ix)}, (14) 
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Table 2: Acceptance probabilities in the "direct" overrelaxation step if the proposed Hip p.'if J4|l 
is w.r.t. the same staple as used in the gauge action {diagonal entries) or another one. These 
figures are specific for the "overlap" projection procedure [T3j as detailed in the appendix. The 
rates have been determined from 100 sweeps on a 8^ lattice, after thermalization. 



with the rescaled staple (normahzation^ such that S^{x) — > 1 if formally the links tend to unity) 
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(15) 



works much better. Typical acceptance rates in this "direct" overrelaxation step at standard 
/?- values (those determined in section 5 to produce matched ensembles) are collected in Table El 
The rates on the diagonal are satisfactory; it is hence advantageous to flip in a quenched SU{3)- 
overrelaxation step with respect to the "effective" staple (|T5|l of the gauge action used. 

The data in this article have been generated with such a mixture of 4 "direct" overrelaxation 
steps per 4-hit Metropolis update, separating configurations by 10 sweeps. While this is prob- 
ably not quite competitive with the standard Cabibbo-Marinari update algorithm |T2], it has 
the great advantage of being easy to implement and represents a reasonable choice whenever 
the CPU-costs are dominated by the measurement. 



3 Smearing details 

The smearing prescriptions considered in this article will be dubbed "APE", "HYP" and 
"EXP" , and they are defined as follows. 

The APE procedure [1^ replaces (in smearing step n) the existing set of links {[7*^"^^^} by 

uirHx) = Psui4 (i-a)f/("-i)(x) + ^s;:^-'\x) } (le) 

where Sjj^~^\x) denotes the 1x1 staple Q built from {U^"'~^^} links. The parameter a governs 
the relative weight of the fluctuation and the projection Psu{3) may be applied optionally. 
The HYP procedure has been introduced in j2] and deflnes the replacement via 

v,M^) = Psui4{i-a,)ui-~'\x) + ^ y: ui-'\x)u;r-'\x+p)ui,''-'\^+f^y} 

U^^\x) = P5^(3){(l-a3)f/("-^)(^) + ^ E KA^)V,A^ + i))K,,{x+fiy} (17) 



^With Psu{3) in dl) as given in the appendix the rescahng is irrelevant, but with an improper one it matters. 
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where I have taken the hberty to interchange ai ^ 03 compared to j2], i.e. our ai denotes the 
fluctuation weight in the first step. Again, ai, as will be considered free parameters and the 
projection will be treated as an option. 

The EXP procedure has been introduced in (there the resulting links are called "stout 
links") and is defined through (no summation over /x) 

[/W(x) = exp (|{[Si"-i)(x)f7i«-^)(x)t-f7^-i)(x)5i"-i)(x)t] - \tT[...]}) Ut'\^) (18) 

where, in principle, any staple could be used, but we shall restrict ourselves to the simplest 
choice, S^j,{x) = S'^^^(x). Upon expanding the exponential, one realizes that this is the re- 
summed version of the ansatz (8) in Ref. The key feature of such a construct is that the 
resulting links t/*^"^ are differentiable w.r.t the input links f/^""^-'; hence the recipe fllHj) may be 
used in a HMC [21111, which is not the case for (fT^ [T7j) unless the projection is abandoned. 



4 Relationship with reduced HQ self energy 

The naive view is that substituting the "thin" links in a Wilson loop with 4D-smeared "fat" 
links replaces the static quark by a quasi-static one which is allowed to rattle inside the 3D-cube 
perpendicular to the HQ line. However, as discussed in ^3], at least for the time-like links in 
the Wilson loop W{R,T) it creates a legitimate - and in fact better - HQ action. 

If the spatial links are subject to 3D smearing, the operator insertion takes place at a specific 
time and the interpretation in the transfer matrix formalism goes through. If the spatial links 
are subject to 4D smearing, positivity may be lost (cf. the short-distance behavior in Fig. 
top), but one can still learn something about the smearing process. If the interpretation [THj 
that smearing reduces the divergent self energy {Sm cxl/a) of the heavy quark is correct, one 
expects that in the ansatz for the expectation value of a Wilson loop at smearing level n 

W^-\R, T) ^ const e-^'"Hi?+T)-<x(")RT (^9) 

the perimeter coefficient gets reduced while the squared string tension stays invariant, i.e. 

^(n) < ^{"-1) (20) 
a^'^) = (21) 

A simple way to test the hypothesis (j^n|) is to look at the T-dependence of an elongated 
(i^-CT) Wilson loop with a fixed R small enough that the perimeter part dominates in (fT^ . 
Such an object may be seen as the 2-point correlator of a "frozen" QQ-state (missing its 
bound-state dynamics), hence a change of its effective mass will only reflect a change of the 
HQ self energies. Fig. |I] (top) shows M^Q{T) = \\og{W{R,T-l) /W{R,T+l)) with flxed i? = 3, 
without smearing and after one step of APE/HYP/EXP smearing, respectively. Obviously, the 
effective mass plateau is lowered [more so under HYP than under EXP or APE] and there is 
a correspondence between this self energy reduction and the factor by which the error-bar is 
decreased and hence the plateau elongated - we shall come back to this point below. 

To test the hypothesis (PT|) it is useful to look at the symmetric diagonal Creutz ratio [12] 
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Figure 1: Top: Mes{T) = llog{W{R,T -1)/W{R,T + 1)) vs. T at R = 3 after one step with 
projection and oape = 0.5,aHYP = (0.3, 0.6, 0.75), q;exp = 0.2. The thick hnks at the initial 
and final time-slice result in a loss of overall positivity. For the relation of plateau-height and 
-length cf. Fig.\^ Bottom: — log(x(a;)) vs. x without smearing and after one step as above. The 
data stem from 500 configurations on a 12^x18 lattice with the Wilson gauge action at (3 = 6.0. 
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Figure 2: Log of noise over signal in the Wilson loop shown in Fig. ^ (top), without smearing 
and after one step of APE/ HYP /EXP -smearing. The full line is a fit to the HYP-data with 
range [10, 14]; the slopes of the dashed lines represent the prediction based on the '^Eqq 
taken from Fig. Q (top) at T = 7. The error on the error is defined as as ((j2— J4)^ + (j4— Je)^)^''^ 
with in the jackknife error at binlength n. Configurations and parameters are those of Fig. ^ 



which (for large x) is independent of fi. In Fig. Q (bottom) — log(x(a;)) vs. x is shown, both for 
the original loop and after one APE/HYP/EXP-step. The squared string tension is invariant 
under smearing, but the error-bar changes. HYP reduces it more aggressively than APE and 
EXP (for the parameters used in Fig. HI the parameter dependence is addressed below). 

In summary. Fig. [T] suggests that the reduced self energy of the quasi-static quark (i.e. its 
divergent piece 5m ~ 1/a) is responsible for the improved signal-to-noise ratio after smearing. 

In fact, this discussion can be made more precise by virtue of an argument due to Lepage [T3] . 
The key idea is rather than determining the variance of a correlator C{x) = (O^(x)O(O)) a 
posteriori from its fluctuation one could estimate it in the MC via (O^(x)O(0)O^(0)O(x)). 
Here, a different symbol has been chosen for the "copy" creation and annihilation operators 
to exclude additional unwanted Wick contractions, i.e. the idea is that they couple to different 
"flavors" (which, apart from the flavor number, have identical properties). In our case this 
noise-squared estimator consists of two Wilson loops, with opposite orientation, on top of each 
other. Asymptotically, it is dominated by two flux-tubes of length zero, i.e. by a state with 
mass 2V{0) = for unitary links. Thus, since the original correlator falls of with an effective 
mass Eqq and its noise is asymptotically constant, the Lepage argument predicts 

signal , , -F -f / N 

— — (t) = const e -^w* 23 
noise 

and Fig. El shows that this expectation is indeed well satisfied. It is clear that the specific 
recipes and parameters in Figs. lH2l are irrelevant - the Lepage argument is a general one. 
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Table 3: Details of the tuning process to get approximately matched lattices. Given is the force 
F{R,T) = l\og{W{R+l,T-l)/W{R+l,T+l)-W{R-l,T+l)/W{R-l,T-l)) on a 8^ lattice 
for R = T = 4. The columns APE/HYP/EXP give the result after one step with the same 
parameters as in Fig. ^ The errors are from a jackknife with binlength 2, using 800 (W) or 200 
(all other actions) configurations per run, after discarding 10 configurations for thermalization. 

5 Matching boxes via the local force 

Our goal is to compare various smearing recipes on lattices generated with either the Wilson 
gauge action or an improved one. In order to give a fair comparison, these lattices should be 
matched, i.e. have the same spacing a in units of a physical scale. Ideally, this is done by 
computing the HQ potential, but rather than moving on to compute various RxT "thin-link" 
Wilson loops, I propose (specifically for the present investigation) two shortcuts. 

The first one is to avoid the T — oo extrapolation and to determine the force F{R, T) = 
l{dR + d*ji)M^g{R,T) with MesiR,T) = ^\og{W{R,T-l)/W{R,T + l)) at one fixed entry, and 
I choose R = T = 4. Since there is no T — oo extrapolation, one can stay in a relatively small 
volume and the finite-volume effect will be the same for all final (tuned) lattices. 

The second shortcut is even less innocent. In order to get precise values for F(4, 4), one 
is tempted to use smeared links. From the observation of Ref. [2j that (one step of) smearing 
affects (for R>2) mainly the additive constant in the HQ potential, one expects F(4, 4) to 
be unchanged, since its arguments are large enough. This is indeed the plan; albeit with the 
proviso that all three smearing recipes will be used to have a cross check. 

Hence, the problem with measuring only F{4, 4) is that the fixed time separation does not 
allow to select a definite state (ideally the groundstate), and the mixture of states that one has 
at T = 4 may vary with the gauge action and the type of link used at T = 0. Therefore, having 
several "thick" links as creation and absorption operators helps, since they couple differently 
to excited states and one may hope that the problem is detected, if it is numerically sizable. 

Table El gives details of the tuning process aimed at finding /5-values which (approximately) 
match a quenched /3w = 6.0 with the Wilson gauge action. The sections "I" and "BKL" refer 
to the tree-level choice in (fTT| IT^ . i.e. ci = ci and C2 = C2 throughout. 
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Table 4: (3-values needed with the Symanzik/Iwasaki/DBW2/Beinlich-Karsch-Laermann/ 
Langfeld actions to have an approximately matched lattice w.r.t. j3 = 6.0 with the Wilson 
gauge action, from interpolating data of Table El The HYP-vaiues will be used subsequently. 

By linearly interpolating, for each action, the simulated couplings one finds the matched 
(according to the local force at R = T = 4 criterion) /?- values. Results are given in Table HI 
without smearing in the first column (where the error is largest) and with one APE/HYP/EXP- 
step in subsequent columns. The unsmeared matched /3-values are compatible with any smeared 
counterpart, but also the agreement among the smeared ones is satisfactory. We rate this as a 
sign that our interpolated matched couplings are essentially correct. Obviously, this does not 
invalidate the theoretical concerns discussed above, but apparently our statistical precision is 
not good enough to really suffer from this. 

In ^\ Necco gives an interpolation formula for log(ro/a) versus (3 for the Iwasaki- and the 
DBW2-action, analogous to the one in for the Wilson gauge action. Solving them, one 
finds that a lattice with /3vv = 6.0 is matched via /5i = 2.635 and /5dbw2 = 1-034 when employing 
one of the RG-motivated actions, respectively. My Pi = 2.652(31) and /5dbw2 = 1.047(14) are 
consistent. In spite of the limitations mentioned above one may thus feel confident that the 
/3-values listed in Table IHyield lattices which are reasonably (say: better than to 10%) matched. 

6 Effect on square Wilson loops and the Creutz ratio 

Whether an improved gauge action is worth the investment is usually decided on the basis of 
the cut-off effects involved. Another practically relevant criterion might be whether such an 
action is amenable to noise- reduction techniques like the APE/HYP/EXP-smearing discussed 
above. Obviously, physical observables should depend as little as possible on the details of the 
smearing process, e.g. the parameters employed. 

Following Ref. [4j, we shall first have a look at a few Wilson loops, comparing the effect of 
changing the parameter (s) and/or the iteration level. 

To this end I have prepared 100 lattices of size 12"^ for each gauge action, using the /3-values 
of the HYP column in Table 01 After this is done, one applies the smearing recipes discussed 
in section 3, varying the parameters in a reasonable range. For HYP I choose to concentrate 
on the direction of the "canonical" a- value introduced in |2j, i.e. I consider six a- values, 

APE : a e {0.20, 0.40, 0.60, 0.80, 1.00, 1.20} 

HYP : a E {0.25, 0.50, 0.75, 1.00, 1.25, 1.50} ■ ag^p with ag^p = (0.3, 0.6, 0.75) (24) 
EXP : a e {0.05, 0.10, 0.20, 0.35, 0.55, 0.80} 

where I deviate from Ref. in my notation aHYP = (astepi, astep2, Osteps)- I choose to apply 2"- 
smoothing steps between successive measurements, resulting in the smearing levels 

ntote{l,3,7,15}. (25) 
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Wilson: c1=0, c2=0; iteration level 1 Symanzik: c1=-1/12, c2=0; iteration level 1 
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Figure 3: For eacii gauge action Wixi, W2X2, ^4x4 plotted versus the parameter number in 
()24|1 ■ indicates the unsmeared starting point. Throughout 1 smearing step (APE, HYP with 
projection) is apphed. Best Wi^i-parameters are oape — 0.7, ohyp — (0.3, 0.6, 0.75), ctEXP — 0.2. 
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Wilson: c1=0, c2=0; iteration level 3 
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Wilson: c1=0, c2=0; iteration level 1 



Wilson: c1=0, c2=0; iteration level 2 
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Figure 5: Wi^i, W2X2, W^4x4 versus the smearing parameter in the APE and HYP procedure 
without projection, using 1 (left) or 3 (right) iterations, on Wilson backgrounds. Note the 
difference to the first graph in Fig.\^ 



Fig.Elshows, for each action, the Wilson loops Wixi, ^2x2, W^4x4 after 1 smearing step (with 
projection back to the gauge group) as a function of the smearing parameter. The plaquette 
gets increased for small a reaching a maximum around a ape — 0.7, ohyp — «hyp; c^exp — 0.2; 
for larger parameters it decreases again. This holds for the larger Wilson loops, too, albeit 
the maximum may be shifted towards larger a (for APE and HYP). What all recipes have in 
common is that excessive parameters drive all loops to (disorder) rather than to 1 (order). 

Fig. m is the same after 7 smearing steps. For small a all Wilson loops are driven more 
efficiently to 1 than before, but beyond a certain "critical" parameter the descent to disorder 
is more pronounced, too. For APE and HYP this culminating a is barely changed compared 
to the previous figure, for EXP it is smaller. 

Fig. El demonstrates a dramatic change in the behavior of any Wilson loop under smearing, 
if the projection in the APE or HYP recipe is abandoned. The graphs have been obtained with 
Wilson glue, but other actions yield the same type of monotonic decrease in a. 

Fig. ini shows the Wilson loops 1^1x1,^^2x2,^4x4 (with projection in the APE and HYP 
cases) versus the iteration level for the second set of parameters in ()24|1 . With such small 
a-values one always observes a benevolent pattern. 

Fig. m shows that with somewhat larger parameters (throughout the fifth one out of the set 
()24|) ) iterating the smearing is not really a good idea. Note that these parameters are to the 
right of the maximum in Fig. El for "sub-critical" a-values no such breakdown has been seen. 

It is worth emphasizing how little the smearing "profile" depends on the gauge action, for 
matched ensembles. A maximum value of P= (Wixi) might be used as a criterion to "optimize" 
the APE/HYP/EXP-parameters (the rather similar values with 1 and 7 iterations are found in 
the captions). Of course, such a criterion is arbitrary, but at least the result is universal, i.e. 
(to a good extent) independent of the gauge action used (for fixed physical spacing). 

More interesting than the behavior of the Wilson loop itself is how the mean and the error 
of a physical observable derived from it evolve under smearing, i.e. whether different recipes 
yield consistent results. For that aim let us consider the squared string tension cr = — log(x) 
defined via the symmetric Creutz ratio (|22I), using the same gauge configurations as before. 

Fig. ISl plots — log(x(5)) versus the parameter a in after 3 smearing steps. Moderate 
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Wilson: c1=0, c2=0; parameter set 2 Symanzik: c1=-1/12, c2=0; parameter set 2 
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Figure 6: For eacii gauge action VTixi, W^2x2, W^4x4 plotted versus the number of smearing 
steps, indicates the unsmeared starting point. Throughout the second parameter in ()24|1 
(APE, HYP with projection ) is used: For small parameters many smearing steps seem favorable. 
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Wilson: c1=0, c2=0; parameter set 5 Symanzik: c1=-1/12, c2=0; parameter set 5 
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Figure 7: For ea.ch. gauge action l^ixi, W^2x2, W^4x4 is plotted versus the number of smearing 
steps, indicates the unsmeared starting point. Throughout the fifth parameter in ()24j) (APE, 
HYP with projection ) is used: For large parameters the result deteriorates after a few steps. 
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Wilson: c1=0, c2=0; iteration level 2 
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Figure 8: For each gauge action the squared string tension a = — log(x(5)) is plotted versus the 
parameter number in ()24|1 . indicates the unsmeared starting point. Throughout 3 smearing 
steps have been apphed. Data with negative central values have been removed. 



15 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



Wilson: c1=0, c2=0; parameter set 2 



APE (proj) 
HYP (proj) 
EXP 

APE (nopr) 
HYP (nopr) 



0.14 
0.12 
0.1 

- 0.08 

- 0.06 

- 0.04 

- 0.02 



Symanzik: c1=-1/12, c2=0; parameter set 2 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



#=1 #=3 #=7 #=15 

Iwasaki: c1 =-0.331, c2=0; parameter set 2 





0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



o APE (proj) 

+ HYP (proj) 

X EXP 

□ APE (nopr) 

HYP (nopr) 



0.14 
0.12 
0.1 

- 0.08 

- 0.06 

- 0.04 

- 0.02 



#=1 #=3 #=7 #=15 

Beinlich-Karsch-Laermann: c1=0, c2=-1/48; parameter set 2 



0.14 



o APE (proj) 

+ HYP (proj) 

X EXP 0.12 

□ APE (nopr) 

HYP (nopr) 



^- 0.1 

- 0.08 

- 0.06 

- 0.04 

- 0.02 



#=1 



#=3 



#=7 



#=15 



o APE (proj) 

+ HYP (proj) 

' EXP 

□ APE (nopr) 

HYP (nopr) 



j=^«-^-:ffla^xfc_- 

:f T---:t__^''<ir - 



#=1 #=3 #=7 #=15 

DBW2: c1 =-1.409, c2=0; parameter set 2 



o APE (proj) 

+ HYP (proj) 

' EXP 

□ APE (nopr) 

<> HYP (nopr) 



-«^^a^^^^ — — — — — — n 



#=1 #=3 #=7 #=15 

Langfeld: c1=0, c2=-0.538; parameter set 2 



o APE (proj) 

+ HYP (proj) 

>^ EXP 

n APE (nopr) 

HYP (nopr) 



—-jr-^M^rM; 1 



#=1 



#=3 



#=7 



#=15 



Figure 9: For each gauge action the squared string tension a = — log(x(5)) is plotted versus the 
number of smearing steps, indicates the unsmeared starting point (where improved actions 
give a smaller error than the Wilson action). Throughout the second parameter in (I24|l is used. 
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parameters suffice to tremendously reduce the noise and and different recipes prove consistent 
within errors. However, as the smearing parameter is raised beyond a reasonable level (typically 
something of the order of the "critical" parameter that achieved the maximum in the plaquette), 
the error may grow again. When repeating this with a higher iteration level basically the same 
behavior emerges, although the onset of systematic distortions with large parameters gets more 
pronounced. There are two more points that deserve a comment. First, I find it astonishing 
that the projection-free APE and HYP recipes manage to decrease the error in a physical 
observable, even though we have just seen that they drive the plaquette to rather than 1. 
Does this indicate that the naive interpretation is wrong, and what one gets is not disorder but 
some type of anti-ferromagnetic ordering ? Obviously, this calls for further investigation, but it 
seems non-trivial to address this issue in a gauge-invariant manner without invoking extreme 
requests for CPU time. Finally, comparing the unsmeared values, one notices that all improved 
actions have an error bar which is smaller than the one with the Wilson action. Hence an 
improved action may lead to a smaller error (with identical statistics) in a gluonic observable, 
but it is clear that this is not competitive with the effect of smearing. 

Fig. iniplots — log(x(5)) versus the number of smearing steps for the second parameter in the 
set (PH). Here, a "plateau" includes 1 to ~7 iterations - further out consistency is lost. With 
a larger parameter a similar picture emerges, but the "plateau" ends earlier. Numerically, a 
value (T~0.06(l) corresponds to a string tension -yo"~490(40) MeV, if a~^ = 2 GeV is assumed. 
Hence a statistical accuracy of 8% can be achieved with 100 lattices of size 12^. 

In summary, these figures indicate a trade-off situation between a statistical error on one 
side and a systematic error on the other. Smearing reduces the former, but when applied in 
excess (either through too large a parameter or too many iterations) it introduces clear artifacts. 
It would be interesting to test how much of these can be attributed to the loss of locality. The 
good news is that a detailed optimization of the parameters is not needed - smearing proves 
beneficial as long as one can explicitly demonstrate that one is in the two-fold plateau; i.e. as 
long as a even larger parameters and/or higher iteration levels do not alter the result. 

7 Effect on field-theoretic topological charge 

Knowing that smearing reduces the UV-fiuctuations and given the relation between UV-fiuctua- 
tions and large renormalization factors and/or cut-off effects, one may ask whether smearing 
reduces the latter. An observable from which one can learn something in this respect without 
having to take the continuum limit is the "naive" field-theoretic topological charge 

qn.i[U] = ^ /^^ tr(Fi2F34 - F13F24 + F14F23) = ^ J dx ti{e^uapF^uF^p) . (26) 

I have implemented ()2fi|) with F^^, defined as the hermitean part of —lUp^p = —il + aF^jy + ... 
where U^y{x) is the mean of the four plaquettes in the (/x, z/)-plane with starting point x. The 
goal is to see what happens if one replaces U^i, by one which is constructed from smeared links. 

Using the same 100 lattices (for each action) of size 12^ as in the previous section, I decided 
to stay with a single parameter fixed at near-optimal values (according to the W{1, 1) criterion), 
ttAPE = 0.5, aHYP = (0.3, 0.6, 0.75), aEXP = 0.2 and three smearing levels ritot ^ {1,3,7}. 

Fig. ^1 shows the time evolution of the bare topological charge (PBjl after 7 smearing steps. 
The APE/HYP varieties without projection tend to drive gnai[f^] to zero for any gauge back- 
ground; these symbols are shown only for the Wilson action. The results of all other recipes 
tend to cluster near non- integer values n/Z with nG Z and Z a function of the smearing recipe 
(the order is 1 < Zhyp < ^ape ^ ^exp for our parameters, i.e. HYP and APE with projection 
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Figure 10: Time history of the bare topological charge ()26p on the matched 12^ lattices after 
7 smearing steps in the definition of F^^ (parameters as in Fig. fTI]) . The projection-free APE/ 
HYP-varieties always give a number very close to zero (omitted in the improved action panels). 
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Wilson, HYP with #=3 and #=7 
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Figure 11: Histograms of the semi-renormalized topological charges after 3 (light) and 7 (dark) 
steps of HYP-smearing. The appropriate Z-factors have been applied, but not the "cast-to- 
integer" operation. The smearing parameter is auYP = (0.3, 0.6, 0.75). 
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Figure 12: Correlation of the unrenormalized topological charges after 1 step of APE- or HYP- 
smearing (light) and of the semi-renormalized charges after 7 steps of APE- or HYP-sniearing 
(dark). In the latter case, and with the final "cast-to-integer" operation performed, the number 
of configurations on which the APE- and HYP-charges disagree is 13, 2, 0, 0, 3, 0, respectively. 
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require only moderate rescaling). This suggests a renormalization factor Z such that 
tends to cluster near integer values. Define - for a given smearing recipe and a given definition 
of Fp,y in ()2(i|l - the charge renormalization factor Z as the solution of 

min(x') where ^ = Y.[z q,,^\\J\- xmvi{Z q^^\lJ\)f . (27) 

Having Z at hand, one may define the renormalized "naive" field-theoretic topological charge 

gren[f/] =nint(Zg,ai[t/]) (28) 

which, by construction, is an integer. Employing such a definition one sees well-defined charge 
histories in each graph of Fig. ^| Obviously, we recover the standard wisdom that - at a fixed 
physical lattice spacing - the Wilson gauge action accomplishes topology changes more easily 
than the Symanzik, Iwasaki and DBW2 action (in this order). The same phenomenon happens 
for the ly^^ ^-actions; the larger C2 the more difficult is it to tunnel into another sector. 

Coming back to the definition ()27| I28 |l . the prime question is whether the assigned charge 
?ren[f^] IS really a property of the gauge background U. Fig. [TTl contains the histograms of Zq 
after 3 and 7 HYP-steps {Z depends both on the smearing level and the gauge action). Here, 
the valleys between the peaks are sufficiently depleted, and the cast-to-integer operation in (PHj) 
makes sense. Unfortunately, there is no real valley structure at the levels and 1. Hence the 
question is to which extent such an integer- valued charge gren[^] is a genuine function of the 
background U and to which extent it is an artifact of the smearing procedure. 

The key to this issue is to compare different smearing recipes. Fig. relates the unrenor- 
malized APE- and HYP-charges after 1 and 7 smearing steps. In the former case there is some 
correlation between the unrenormalized gnai[^], but there is no clustering whatsoever. In the 
latter case, i.e. after a few more smearing steps, the two rescaled Zq^^i cluster predominantly 
near the same integers. With 100 configurations in each ensemble, on 87/98/100/100/97/100 
the 7-APE-step and 7-HYP-step charges agree after renormalization. Thus, for ~ 2 GeV 
and an improved gauge action the two charges agree almost on a configuration-by-configuration 
basis and this suggests that such a charge is indeed a function of the background U. 

Still, this does not mean that the renormalized "naive" field-theoretic charge ^I7\ 0H|) is as 
much a clean concept as the fermionic charge defined through the index of the Dirac operator 
Qind = |tr(75 D) [ini EO] with D any Dirac operator that satisfies the Ginsparg- Wilson rela- 
tion [21]. However, ()27l I28|) has the virtue of being cheap in terms of CPU-time; the 0(3. ..7) 
smearing steps cost less than the hundreds of cooling steps that have been used in the past, and 
results are competitive with those of the cooling era (see e.g. Ref. [22] for a nice example). An 
issue that deserves further study is whether too much smearing eventually destroys the charge. 

8 Summary 

The goal of this work has been a comparison of the newly proposed EXP smearing ^ against 
the well established APE and HYP procedures 12] on various gauge backgrounds. 

I started by scanning through (part of) the space of parameters and iteration levels to see 
which combinations drive the plaquette and larger Wilson loops to 1 rather than and to see 
whether the choice of the gauge action would have a big infiuence. The result is that standard 
parameters around oape — 0.65, ohyp — (0.3, 0.6, 0.75), oexp — 0.15 prove nearly optimal, at 
least if 0(1. .3) iterations are performed. Obviously, the plaquette criterion is rather arbitrary, 
but for a^^^2 GeV the result is universal, i.e. independent of the gauge action. 
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The only physics point of this paper has been an exphcit test of the Lepage argument in 
the heavy-heavy case. This argument estabhshes a connection between the signal-to-noise ratio 
of an observable involving a static quark and the HQ self energy, and our Fig.|21 shows that it 
passed the test with flying colors. Hence a future optimization strategy tailored to HQ physics 
might be to reduce the HQ self energy as much as possible. The key insight is that in the case 
of large Wilson loops various smearing recipes differ in the signal generated, not in the noise. 

Next I evaluated the squared string tension defined via the Creutz ratio at i? = T~ 0.5 fm. 
The main finding is that any smearing recipe drastically reduces the noise and a detailed tuning 
of parameters and/or iteration level is neither needed nor useful. Instead there is a broad two- 
dimensional plateau in which changing the two does not significantly alter the central value or 
the error. Only at absurdly high levels systematic effects get visible, presumably due to the 
reduced locality in terms of the original links. In summary, moderate smearing suggests itself 
as a simple and efficient means to damp the UV-noise in standard gluonic observables. 

Finally, I have focused on the "naive" field-theoretic topological charge. Here, the power of 
any smearing recipe depends a lot on the gauge action used. The good news is that the lattice 
spacings where the smearing details implicit in the definition ()261 071 OHj) prove irrelevant are 
accessible; for the Iwasaki action it is around a~0.1fm, for other actions not much below. 
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well as the benefits from a discussion with Rainer Sommer, Michele Delia Morte and Francesco 
Knechtli on the Lepage argument. Computations were done on a PC at DESY Zeuthen. This 
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Appendix: Projection details 

Here I spell out the projection procedure introduced in P^. For related ideas see j^HlEl!- 
The idea is to define the projection step V^(x) Wf^{x) via the overlap-type recipe 

' [det (numerator)] V3 ^^^^ 

which leads to an eigenvector decomposition problem of the hermitean operator X = V^V. 
In principle, one might find the eigenvalues xi, ...,X3 of X (following Cardano) by forming 

a = -tr(X) 

b = X11X22 + X22X33 + X33X11 — X12X21 — X23X32 — X31X13 
c = -det(X) 

V = b-j 

2a^ ab 

and continuing with 

r if g = ^ xi = Ml - f 

ifp = Othenui = J (l + i^^)^) ifg>0 [and X2 = (-^ + i^)ui - f (30) 

ifg<0j a;3 = (-i-if )Mi-f 



3/ / 2 3 Xi = Ml + 171 — 3 

else = V~2 + V^ + It and X2 = (-^ + i^)ni + (-i - i^)Mi - f (31) 
^1 = ^3 = (-i - if )wi + (-i + if )^i - f 
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where the roots in (j3Tjl must be taken over the body of complex numbers, though the result 
xi, X3 is real. 

A practical problem emerges whenever X is close to unity (what is generic after a few 
smearing iterations), as the naive implementation cannot separate the eigenvalues any more. 
In order to reduce numerical extinction, one should implement 

Y = ^{V^-l){V-l)+h.c. + iV^-l) + iV-l) { = X-l) (32) 

rather than X and trade, for the same reason, the quantities a, b, c for alternative ones designed 
to involve only "small" numbers in the case X ~ 1 : 



as = -Yu-Y22-Ys3 ( = a + 3^0j (33) 
63 = YnY22 + Y22Y,; + Y,;Y^,-YuY2i-Y23Y,2-Ys,Yrs-2as ( = 6-3 ~ o) (34) 

Cl = —^11^22 ^33 ^ ^12^23X31 — ^32^21^13 

+YnY2A + Y22YsiY,s + ^33^12^21 - h - a-, ( = c+ 1 ^ o) (35) 

P = &3-y + 2a3 (36) 

q = ^-^ + «3-^ + &3 + ci (37) 

before one proceeds with (jHUll^ . where again yi = ui — y, etc. should be formed rather 

than xi,X2,X3. Knowing the eigenvalues (minus 1) with good precision, it is easy to find the 
eigenvectors of X by solving the characteristic equation 

Yv^ = yiV, («G{1,..,3}) (38) 

and after normalizing them one forms 

W = . (39) 

[det (numerator ) J ^/ 

Of course this costly projection procedure is only needed if the argument is potentially far 
away from SU{3); whenever V is known to be close to SU{3) (e.g. in the MC a newly proposed 
link is projected to SU{3) before it is subject to the accept /reject-test, but here the deviation 
is just due to round-off errors) one may project via an approximate procedure, e.g. via 

V V' = lv- \VVW , V V" = V'(l- \{dei V - 1)) . (40) 

Finally, it is worth mentioning that the same eigenmode decomposition technique (simplified 
by a = 0) has been used in the matrix exponentiation in the EXP recipe ()18|1 . 
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